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ABSTRACT 

The Crab pulsar has suffered in 1975 and 1989 two glitches in which the frequency 
did not relaxed to the extrapolated pre-glitch value but rather spun up showing long-term 
changes in the frequency derivative O. This particular behaviour has been interpreted 
as evidence for an evolution of the torque acting upon the star. A variable torque may 
be related to non-canonical braking indexes, for which some determinations have been 
possible. We briefly analyse in this work the consistency of postulating a growth in the 
angle between the magnetic moment and the rotation axis as the cause of such events. We 
show that this hypothesis leads to the determination of the initial period, initial and present 
angles, according to the assumed angle growth, for young pulsars whose respective braking 
indices Uobs and jerk parameters mobs are known, and some insights on the equation of 
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state. 
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1 INTRODUCTION 


The dynamics of compact stars offers a unique opportunity to explore and understand 
the physics of dense matter. Particularly, timing irregularities have proved to be extremely 
important (and challenging) for building a coherent picture of those stellar interiors 
involving hadronic matter at sub nuclear and supranuclear densities. 

Among the most notorious timing irregularities, sudden discontinuities (glitches) are 
observed in the pulsar frequency O and spin-down rate O. Several of such events have been 
observed from the Crab, Vela and some other pulsars. It is generally agreed that they are 
the result of a (complex) interplay between the superfluid component(s) and the rest of the 
star. Theoretical models searched an explanation in terms of starquakes (Baym et al.l969) 
or, more recently, vortex motion in the superfluid component (Alpar et al. 1984, Pines 
& Alpar 1985, Link, Epstein & Baym 1993). While the latter model is flexible enough to 
accomodate a large body of observations, evidence coming from very timely and detailed 
observations of two glitches in the Crab pulsar seems to indicate that new physical inputs 
may be required to explain the data (see below). Since the Crab can be considered as the 
best studied object, it may be presumed that an analogous behaviour could be present 
in the young pulsar sample for which an increasing body of data is available. We shall 
attempt to understand the dynamics of the youngest objects and build a consistent picture 
of their evolution. 

In 1975 and 1989 the spin rate of the Crab pulsar O suddenly increased by amounts 
AO/fl ~ 10“® and after that continued to spin-down at a faster rate (that is, the pulsar 
continued to spin slower than before) of AO/0 10 ^ (Gullahorn et al. 1977, Lohsen 
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1981, Lyne, Smith & Pritchard 1992). The same feature is also present in the 1969, 1981 
and 1986 events (Lyne & Pritchard 1987, Lyne, Pritchard & Smith 1993). This behaviour 
has been attributed to a decoupling of some internal shell or a change in the external 
torque acting on the star (Gullahorn el al. 1977, Demiahski & Proszyhski 1983). By using 
a general postglitch relaxation equation of the form 


du{r, t) 
dt 




for the differential rotation w as a function of the specihc external torque f{u),r), Link, 
Epstein and Baym (1992) have argued that a frequency dehcit can not be explained by the 
existing glitch models. In other words, according to their work the postglitch frequency 
flc(t) must be always greater than the extrapolated preglitch frequency Oco(^), contrary to 
the observations. On the other hand, Alpar & Pines (1993) have discussed the theoretical 
interpretation of these events in the framework of vortex creep theory. Even though there 
seems to be enough room for such a behaviour in the latter, several details of the pinning 
layers are still unclear and complicate the interpretation. Thus, it may be interesting to 
explore alternatives to bring the theoretical picture closer to the observed phenomenology. 

An important aspect of these observations is that, given the much larger amount 
of data taken from the Crab pulsar and the difficulties of extracting this signature from 
timing noise, it is entirely possible that other young pulsars also behave similarly. If so, 
events producing a permanent AO/0 may be important for the understanding of pulsar 
dynamics. Models having variable external torques have been considered in the past (see 
Blandford & Romani 1988 and references therein) and may be helpful for understanding 
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the data. Our goal in this work is to see to what extent a simple (but consistent) picture 
of young pulsar dynamics can be built by assuming a specihc version of a variable external 
torque model. Section 2 is dedicated to formulate a minimal dynamical model in which 
the angle between the magnetic dipole and the rotation axis of the pulsar is allowed to 
vary according to simple laws. Section 3 presents an application of the model to the Crab, 
Vela and two other interesting pulsars (PSR B1509-58 and PSR B0540-69 ). Finally, we 
present a discussion and conclusions in Section 4. There is also an Appendix with exact 
solutions for the proposed growth laws. 

2 BASIC EQUATIONS 

The equation of motion of a rotating pulsar, assumed to have a crust+core normal 
component (suffix c) and a superfluid shell (suffix s) is 

Ic^c T Is^s — '^ext (1) 

since the frequency difference Oc ~ bis = is constant on sufficiently long timescales, the 
rotational equilibrium implies further that 


Itot^c — “Text- ( 2 ) 

Up to now studies have been based mainly on the use of torque laws of the form 
Text = — where A is a constant (up to its possible decay for very old objects) 

and the braking index n is assumed to be 3 as predicted by the magnetic dipole model 
(see e.g. Manchester & Taylor 1977). The vacuum dipole model expression corresponds 
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to K = \M\^ sin^ a, where a is the angle between M and , |M| = So and Bo 

, R are the magnetic held and the radins of the star respectively. 

From eq.(2) and the form of Text it is clear that the pecnliar events of 1975 and 1989 
reqnire either a redaction of Rot, an increase of the magnitnde of the magnetic moment 
M or an increase of the angle a between M and flc- Althongh all these hypothesis are in 
principle allowed by the data, we shall address here the change in the relative orientation 
of M and flc axis in those events. The idea is appealing not only becanse of its simplicity 
bnt also because it may be considered as a realisation of Ruderman’s plate tectonic theory 
(Ruderman 1991) where it hnds a natural room. The other two hypothesis also have a 
theoretical support in the framework of vortex creep theory (Alpar & Pines 1993) and 
magnetic held generation or surfacing (Blandford, Applegate & Hernquist 1883, Muslimov 
& Page 1996) respectively. They have been previously considered and we shall not address 
them here. 

At this point it must be noted that a counter-aligning pulsar in which a increases 
with time is quite remarkable, since the opposite situation is to be expected naively for 
the rotating/radiating star. Therefore, the presence of ’’anomalous” glitch events is per se 
a signature of the complexity of the internal/magnetospheric dynamics which ultimately 
determine the behaviour of a{t). Several papers have addressed the question of the internal 
dynamics and its consequences for a{t) (Michel & Goldwire 1970, Davis & Goldstein 
1970, Goldreich 1970, Michel 1973, Lyne & Manchester 1988, Michel 1991). Macy (1974) 
has given a general treatment by solving simultaneously the system of equations which 
determine a{t) and Dc(t) (= ^puise) adopting a simple vacuum dipole model for the pulsar 
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radiation. In order to check the consistency of the increasing angle hypothesis we have 
chosen instead to parametrise the growth of a by simple expressions and solve analytical 
models to compare them with the observations. 

In hrst place, we want to know if the angle growth happens solely on glitches, in a 
discrete mode, or if there is a inter-glitch continnons contribntion. To do this, we can 
determine the variation of the angle a dne to persistent shifts AO/0 and AO/0 which is 
easily fonnd from eq.(2) to be 


Act 


AO ^AO^ tana 

H' ~ ) 2 


(3) 


Dividing eq.(3) by At, here dehned as a typical time-scale between glitches, we obtain 
a (mean) increase rate which we shall denote as (^). We shall keep in mind that, since 
all we have at disposal is a very short observational span of ~ 20 yr at most, there conld 
be a discrepancy between (^) and the continnons modeling in which d is a trne local 
derivative. In fact, the data from Lyne, Pritchard & Smith (1993) allows ns to hnd for the 
Crab pnlsar 


/ Aa—1—\ ^ i_8 X 10 ^ radyr ^ (4) 

\ At tana/ ^ ^ 

with At = 4.6 yr. On the other hand, the rate needed to acco un t for the observed braking 

index is fonnd from eq.(Al4) to be 


a 


tana 


'IT'obs 3 O 
2 O 


= 9.6 X 10 ^ radyr 


-1 


(5) 


abont 5 times the rate obtained from glitches. For PSR B0540-69, Vela and PSR B1509-58, 
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the rates given by eq.(5) are respectively 14.4 x 10 ^ radyr 3.49 x 10 ^ radyr ^ and 
2.62 X 10 ^ radyr All these valnes are within one order of magnitnde, strengthening 
the idea that they conld have the same origin. Once we have acknowledged that the main 
contribntion to the angle growth comes of the interglitch slowdown, and becanse even 
pnlsars (like PSR B1509-58 and PSR B0540-69) which did not display any glitches nntil 
now show low braking indexes, we have chosen to describe the angle growth as a continuous 
fnnction rather than a discrete one. The other hypothesis pnt forward to explain the low 
braking index conld be invoked here as above, bnt we will only explore the angle growth 
scenario. 

We have tried fonr simple cases for the angle growth, namely an exponential, linear, 
power-law and logarithmic fnnctions (see Appendix). Analytical expressions for 0(t) have 
been obtained in all cases which, in principle, may be acceptable as descriptions of the 
long-term pnlsar’s spindown. However, a closer inspection of these fl{t) forms reveals that 
the linear and power-law solntions are not of general applicability becanse they critically 
depend on the specihc inpnt parameters (like the present angle ctp) to exist and be positive 
dehnite. Therefore, we have been led to consider the exponential growth a{t) = cto 
and the logarithmic growth a(t) = In [ A (^g«p _ _l_ g«o j -which do not snffer from 

these misbehavionrs (see Appendix). 

We now tnrn to the application of a varying-angle model described by the formnlae 
of the Appendix to specihc cases in the next Section. 



3 APPLICATION OF THE MODEL 


3.1 General Considerations 

In order to check the consistency of the varying-angle model we shall apply it to 
the fonr pnlsars in which O (and hence Uots) has been measnred. Starting from the Crab 
pnlsar, in which a rather direct measnre of mobs has been also obtained, we have attempted 
to nnify as far as possible the behavionr of the yonng pnlsar sample. The strategy is as 
follows: assnming that the braking index and jerk parameter discrepancies are solely dne 
to the growth of a, we impose the observed valnes of nobs and mobs (see eqs.(Al7) and 
(A19) in the Appendix) to determine a solntion for a. Once these valnes are obtained, 
eq.(Al), (A2) and (A3) can be nsed to calculate Po, cto and the time-scale for maximum 
torque (that is, the age at which a reaches 7r/2) termed tm- 

3.2 The Crab Pulsar (PSR B0531+21) 

The Crab pulsar has been extensively studied and the basic quantities like nobs and 
mobs determined through measurements of O derivatives. Using the observational data at 
MJD 40000.0 available from Lyne, Pritchard & Smith (1993), we found for the exponential 
form (the logarithmic form does not provide a consistent solution for this pulsar) 

ap ~ 68 ° 

tto ^ 56° 

Po — 19 ms (6) 
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tm - 2300 yr. 

We note that the tm obtained is larger than the trne age of the Crab tp ~ 940 yr, as 
it shonld be. It is also worth mentioning that ctp is more consistent with the new valne 
inferred from optical polarization data (ctp ~ 60°, F.G.Smith et al. 1988) than the older 
one (ttp > 80°, Kristian et al. 1970) and radio data (Rankin 1990). 

As a corollary of the dynamical solntion, the strnctnral constant C — 2 ,c^°it t 
Appendix) can be now evalnated to be 4.14 x 10“^® s. If the magnetic held is taken 
to be ~ 3.8 X 10^^ G (Taylor, Manchester and Lyne 1993), a constraint on the prodnct 
(( 1045 gcm 2 ) ^ — 1.16is obtained from the torqne expression. We thns 
see that the proposed dynamics offers some insight onto the eqnation of state, namely, 
the relationship favonrs a eqnation of state slightly softer than the Bethe-Johnson I model 
(Bethe & Johnson 1974). 

3.3 The Vela Pulsar (PSR B0833-45) 

A measnrement of Uobs for the Vela pnlsar has been recently obtained by Lyne et al. 
(1996), and it has not been possible determine mobs as yet, therefore the procednre nsed 
above is not snitable. However, introdncing the observational data from the catalog of 
Taylor, Manchester & Lyne (1993), we hnd that a logarithmic solntion is possible only if 
mobs > 4.5 (and O > —8.6 x 10“^^ rads"^), and an exponential one reqnires mobs < 3.2 
(and Q < —6.1 x 10“^^rads“^). An observational determination of the jerk parameter 
wonld thns discriminate the angle-growth. Meanwhile, eq.(Al5) provides an npper limit 
for mobs of 3.2, calcnlated without considering the term which contains the angle and its 
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derivatives. Therefore the only acceptable solution is a exponential growth of a. For an 
arbitrary value of ctp we have found the limits 

cto < 32° 

Po > 51ms (7) 

tm < 1-4 X 10^ yr 

where the carachteristic age has been used as the true age. There are claims (Aschenbach, 
Egger & Triimper 1995, Lyne et al. 1996) suggesting that Vela is actually older by a factor 
of 2 or 3. This last possibility would imply in our analysis smaller limits for cto and Po for 
the same ctp, and a bigger limit for tm ■ An older pulsar also implies that mobs becomes a 
positive number. This is as far as we can go without any further data. 

3.4 PSR BO54 0-69 

This pulsar does not have a observational determination of mobs either. As an 
interesting feature of our models we predict using the data presented in Taylor et al. 
(1995) (Taylor, Manchester & Lyne 1993) a logarithmic angle growth if mobs > 7.8 (O 
> —8.4 X10“^^ rad s“^), and an exponential one if mobs < 7.2(0 < —7.7x10“^^ rad s“^). 
The upper limit from eq.(Al5) is mobs — 7.3, so a logarithmic growth is also excluded as 
in the case of Vela. Letting ctp to be a free parameter, we can calculate for an exponential 
growth 
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cto < 41° 


Po > 23 ms. (8) 

tm < 3.3 X 10^ yr 

3.5 PSR B1509-58 

With the data by Kaspi et al. (1994) and using the characteristic age as tp ~ 1550 yr, 
we hnd that the logarithmic growth is the only one allowed for this pulsar and gives 

ap ~ 83° 

cto ~ 61° (9) 

Po — 45 ms 

~ 2200 yr. 

Again, we can calculate the structural constant (7 ~ 6 x 10“^^ s, which leads to 
i iQ^cm T ( io4&gcmO ^ “ 1-01 S'assumiug Bo ~ 1.55 x 10^^ G (Taylor, Manchester 
and Lyne 1993). This is approximately the same result obtained for the Crab pulsar. 
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or about 3% relative difference in the radius of both objects, suggesting that the same 
equation of state can be applied to them. 

There is some dispute about this pulsar’s true age (Thorsett 1992, Blandford & 
Romani 1988, Van den Bergh & Kramper 1984). If our picture is correct, the maximum 
admissible value for the pulsar’s age tp is ~ 1730 yr but in this case the initial period 
would have been as low as ~ 4 ms, which seems difficult to accept. These hgures would 
support the association of this pulsar with the 185 AD supernovae. On the other hand, if 
this pulsar is older (the remnant is calculated to be 10^ yr, see Van den Bergh & Kramper 
and references therein), its angle a may have reached 7r/2 in the past, so that our analisis 
no longer applies to it. 


4 CONCLUSIONS 

We have studied simple parametric dynamical models of the angle growth which led 
us to determine the initial features of the four considered pulsars. Even if we have tried 
to unify as much as possible their dynamical behaviour, some ambiguity related to the 
intrinsic difficulty of the data analysis remains. As an example, the data analysis of the 
Crab pulsar made by Lyne, Pritchard & Smith (1993) have not directly calculated the 
third derivative of frequency u , but instead proceeded to End it by setting mobs = ^ 
j = nobs{‘^'nobs — 1). A inspection of eq.(Al5) reveals that if the canonical braking 
index is replaced by nobs (M.P.Allen & J.E.Horvath, in preparation) the angle-dependent 
term automatically vanishes, but since the variation of Uobs is small in the Crab, the 
correction is also small. Strictly speaking, the true mobs in the case of a varying angle 
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would be greater than that value if the ” angular” term 


Vt f a 2a \ , , 

7^\~ ~ • /o N < '>^obs - 1 (10) 

O \a sm( 2 a)/ 

and smaller if the inequality is reversed. If we substitute ctp = 68 ® as found in Section 
3.2 in eq.(Al5), we obtain a correction to mobs of + 0.15. It should be noted that this 
correction in O implies small corrections in the other derivatives measurements particularly 
in O. Also, the discrepancy in the braking index, either originated by angle growth, 
magnetic held increase or inertia moment decrease, necessarily leads to some correction 
in the estimatives of mobs rnade through the simple assumption above, unless fortuitous 
cancellations occur. 

We consider that an exponential angle growth, with a typical (but non-universal) time- 
scale < 10^ yr is a good (yet simple) model for the young pulsars dynamics. It accounts for 
the observed discrepancies in the braking index and the jerk parameter. It also provides 
limits for the initial parameters of individual pulsars if just Uobs is known. These preditions 
can be compared to the observational data to be gathered in a few years. 

As a hnal remark we stress that our results show that the characteristic age may be 
a rather poor estimative of the true age. The dependence of the torque on t forces a full 
calculation of tp by inverting the complicated expression of eq.(Al) instead of using the 
standard form 








n—1 


( 11 ) 


(n — l)Op 

which nevertheless remains valid (even if d 7 ^ 0) in Uobs = constant (M.P.Allen & 
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J.E.Horvath, in preparation). We conclnde that, even if a exponential growth of the a 
angle seems to be consistent with several observed aspects of yonng pnlsars behavionr, 
more observational data (mainly a detection of a change in the average pnlse prohle at the 
level ~ 10“^) is needed to formnlate more complete and complex models of the dynamics. 
Reliable determinations of a for these pnlsars throngh optical polarization or other methods 
and niobs wonld be decisive to test angle-growing models. 

While this paper was being written we have received a preprint by B.Link & 
R.I.Epstein where some of these ideas were independently exposed and discnssed. 
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7 APPENDIX 


We present in this Appendix the explicit forms a{t) employed for the determination 
of the dynamical pnlsar histories. First let ns consider an exponential growth of the angle 
a{t) = aoe*/*“. Assnming a vacnnm dipole model for Text^ eq.(2) can be integrated to 
yield the analytic form for O 


0 = 0 , 


where T = — ^ = ((7 sin^ ^ is a slowdown time-scale at the present time, with 

C being a strnctnral constant depending on the chosen eqnation of state, ctp and cto are 
the present and initial angle respectively, and = apjap is a constant time-scale for the 
growth of the angle. Now, we can express the time-scale for maximnn torqne the initial 
period Po and cto in terms of observable qnantities yielding 


1 + 


tr 




{2apy^ - (2a„) 


2n n -1/2 


Tsin^ap^^^ (2n)(2n)! 


, -1 


(Al) 


t 


m 



+ tp 


(A2) 


— m tp/io 


OLo = OLpt 


(A3) 


If we assnme instead a linear growth law a{t) = (§ — cto) y—^ and integrate again 

the eqnation of motion, we obtain the non-linear algebraic eqnation 


0 = 0 , 


1 - 


1 


Ta sin^ a. 


ap — a — sin(2ap) -f sin(2Q;) 


- 1/2 


(A4) 


which gives 
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t 


m — 




1 

— tp 
a 


{Ah) 


oio = oip — atp. 


(-46) 


Integration of a logarithmic law a{t) = In ^ — e“°) + e“° gives 


O = np<' 1 - 


2tr 


hT{e^p — e“°) sin^ 


X 


X 


e«P 


sin^ ar 


— sin( 2 ap) + 2 ) — e 


sin' 


a — sin( 2 a) + 2 


- 1/2 


(A7) 


with the following resnlts 




— ctp — In (l — tpCXp) 


e“p Ctr 


(A8) 


olq — rip “ 1 “ In tp hp^. 


(A9) 


Finally, integrating a power-law a{t) = (f — + cto , where 6 is a positive 

constant 7 ^ 1 , we have 


O 



1 - 


tp tm 


rri / TT 1 • 2 

J- [2 — ctoj Sin a 




! X 


X 


cos(2ap) E /2i+2(«p) 

i=Q 


cos( 2 a) E /2i+2(a) + 

i=Q 
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2b 2 ) 

+ sin(2Q;p) E f2i+l{ 


Ctr 


sin(2a) ^ / 2 i+i( 


a] 


j=0 


i=0 


where 


-1/2 


(AlO) 


/ 2 i+i(a) — (- 1 )* 


{a-ao)^ ^ 

22i+i(i _ 2i - 1)! 


/ 2 i+ 2 (a) = (- 1 ) 


i («-««)'> 


T-2i-2 


22*+2(i - 2i - 2)! ’ 

and int{x) is the greatest integer nnmber lesser than x, for h = -I, ■■■}■ For any 

other b we will have 


n = 1 


tp tm 


rj-i / TT N • 2 

J- [2 ~ cto) Sin a 


(ctp — cto) (cos^ cto — cos 2 ap) — 


[a - 


cos^ Ctr 


cos 


2 a) 



00 

cos(2ap)^5f2i(ap) 

i=l 


00 

cos( 2 a)^ 5 f 2 i(a) + 
i=l 


00 00 

+ sin(2ap)^5f2i+i(ap) - sin(2a)^5f2i+i(a) 
i=0 i=0 



(All) 


where 


j 22 *(a - ao)b+ 2 * 

»«(“) = (-1) a + a)! 


^ 2 i+i(a) = (- 1 ) 


;22"+i(a - ao)^+2"+^ 

(5 + 2; + !)! 
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giving 


t — 

Lm — 


(f -ap)6 


, 1/6 


Ctr 


(- 412 ) 


(Xq Xp 


tpOip 


(-413) 


Note that the expressions eqs.(Al), (A4), (A7), (AlO) and (All) mnst be now nsed 
to dehne the characteristic age of the pnlsar replacing the nsnal form r = , t^t (see 
Section 4). 

It also follows that in these angle varying models it can be shown that, even if the 
power of O in the torqne expression is exactly 3, the observed braking index dehned as 
Uobs = picks np an extra term 


^obs 3 2 


n 


da 

cot a —- 
dt 


and the observed jerk parameter dehned as mobs 


Q is also modihed 


(- 414 ) 


miobs 


3 (2 TT-obs 1) A ij^obs 3) 


o 

^obs A '77 



2 d \ 
sin(2Q;) J 


Therefore, nsing the simple expressions for a{t) we obtain 


(- 415 ) 


sin(2Q() 2 q( TTlobs f^obsij^obs T 3) T 3 
4 cosset {nobs — ‘^Y 


(A16) 


for the exponential form. 


cos a 



miobs '^obsi^'IT'obs T 3) T 3 
{nobs - 3)2 


- 1/2 


(A17) 
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for the linear form, 


sin(2(r) + 2 _ mobs - ngbsiriobs + 3) + 3 
4 cosset {riobs — ^Y 

for the logarithmic form, and 


(A18) 


cos a 


2 I '^obs '^obsijf'obs 4“ 3) -|- 3 ^ T (h 1) 

ij^obs 3)^ tp {jlobs 3) 


- 1/2 


(A19) 


for the power-law. These forms express the nnknown angles in terms of the Uobs and mobs- 
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